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Abstract 

We prove the all-order exponentiation of soft logarithmic corrections to prompt photon pro- 
duction in hadronic collisions, by generalizing an approach previously developed in the context of 
Drell-Yan production and deep-inelastic scattering. We show that all large logs in the soft limit 
can be expressed in terms of two dimensionful variables, and we use the renormalization group 
to resum them. The resummed results that we obtain are more general though less predictive 
than those proposed by other groups, in that they can accommodate for violations of Sudakov 
factorization. 
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1 Introduction 



The resummation of logarithmically-enhanced terms in the perturbative expansion of QCD cross 
sections near the boundary of phase space (Sudakov [T] resummation, henceforth) has been 
proved to next-to-leading logarithmic order for deep-inelastic scattering in the x — > 1 limit and 
for Drell-Yan production near threshold [21 Ej. In fact, resummation formulae can be derived for 
a wider class of processes and to all logarithmic orders, by assuming the validity of a suitable 
two-scale factorization 0j. The two scales correspond to the hard scale which guarantees the 
perturbative nature of the process, and another dimensionful variable, whose large ratio to the 
first perturbative scale is the quantity whose log must be resummed to all orders. This two-scale 
factorization, originally derived for semi-inclusive two-jet production in e + e~ annihilation [5], 
can be proved for deep-inelastic scattering and Drell-Yan in the 3 theory in six space-time 
dimensions |6J, but its generalization to QCD processes is nontrivial. 

In ref. [7j a new approach to Sudakov resummation was developed, and applied to the Drell- 
Yan and deep-inelastic scattering processes. This approach has the advantage of being valid to 
all logarithmic orders, and self-contained, in that it does not require any factorization beyond the 
standard factorization of collinear singularities. It relies on an essentially kinematical analysis of 
the phase space for the given process in the soft limit, which is used to establish the result that 
the dependence on the resummation variable only appears through a given fixed dimensionful 
combination. This provides a second dimensionful variable, along with the hard scale of the 
process, which can be resummed using standard renormalization group techniques. Beyond the 
leading log level, the resummed result found within this approach turns out to be somewhat less 
predictive than the result of refs. [21 E] : m those references the next fc -to-leading log resummed 
result is fully determined by a fixed next fc -to-leading order computation, whereas a higher fixed 
order computation is needed to determine all coefficients in the resummed formula of ref. [Zj. 
The more predictive result is recovered within this approach if the dependence of the perturba- 
tive coefficients on the two dimensionful variables factorizes, i.e. if the two-scale factorization 
mentioned above holds. 

In this paper we discuss the generalization of the approach of ref. |7J to the resummation of 
the inclusive transverse momentum spectrum of prompt photons produced in hadronic collisions 
in the region where the transverse momentum is close to its maximal value. Prompt photon pro- 
duction is a less inclusive process than Drell-Yan or deep-inelastic scattering, and it is especially 
interesting from the point of view of the approach of ref. [7], because the large logs which must 
be resummed turn out to depend on two independent dimensionful variables, on top of the hard 
scale of the process: hence, prompt photon production is characterized by three scales. The 
possibility that the general factorization ref. [U] might extend to prompt photon production was 
discussed in ref. jS] , based on previous generalizations 9 J of factorization, and used to derive the 
corresponding resummed results. Resummation formulae for prompt photon production in the 
approach of ref. [2] were also proposed in ref. JUj > and some arguments which might support such 
resummation were presented in ref. ^T] . Our treatment will provide a full proof of resummation 
to all logarithmic orders. Our resummation formula does not require the factorization proposed 
in refs. [El EI], and it is accordingly less predictive. Because of the presence of two scales, it 
is also weaker than the result of ref. [7j for DIS and Drell-Yan production. Increasingly more 
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predictive results are recovered if increasingly restrictive forms of factorization hold. 



2 Kinematics of prompt photon production in the soft 
limit 

We consider the process 

H 1 (P 1 )+H 3 (P 2 )^ 7 (p 1 )+X 1 (2.1) 

specifically the differential cross section p 3 ± -^-, where p± is the transverse momentum of the 
photon with respect to the direction of the colliding hadrons Hi and H 2 , and 

* = ^p S = (P 1 + P 2 ) 2 ; 0<x ± <l. (2.2) 

The factorized expression for this cross section in perturbative QCD is 

da f 1 ( Q 2 \ 

p\- — {x,p\) = V / dx 1 dx 2 dzx 1 F^ 1 (x 1 ,/j 2 ) x 2 F" 2 (x 2 , /j 2 ) C ab z, — , a s (/i 2 ) \8{x- zxix 2 ) , 
dp± ^Jo V V- ) 

(2.3) 

where F^ l {xx, ji 2 ), F b H2 (x 2 , /x 2 ) are the distribution functions of partons a, b in the colliding 
hadrons, we have defined 

Q 2 = ±P 2 ±, (2.4) 
O 2 

z = — , 0<z<l, (2.5) 
s 

where s is the center-of-mass energy of the partonic process, and the coefficient function 
C a b(z, %,« s (a* 2 )) is defined in terms of the partonic cross section for the process where par- 
tons a, o are incoming as 

=,*fS. (2.6) 

We will study the cross section eq. ()2.3|) in the limit in which z — > 1, i.e., the transverse 
momentum of the photon is close to its maximal value, thereby suppressing the phase space for 
further parton radiation. The convolution in eq. ()2.3j) is turned into an ordinary product by 
Mellin transformation: 

a(N, Q 2 ) = £ (N + 1, /x 2 ) F b H \N + 1, /i 2 ) C ab (N, Q 2 /fi 2 , a s (/i 2 )), (2.7) 

a.b 

where 

a(N,Q 2 )= f X dxx^pl^-ixipl) (2.8) 

JO Ctp_|_ 

and similarly for C a b and Ff ] . 

Whereas the cross section a(N, Q 2 ) is clearly /^-independent, this is not the case for each 
contribution to it from individual parton subprocesses. However, the fi 2 dependence of each 
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contribution to the sum over a, b in eq. (|2.3|) is proportional to the off-diagonal anomalous di- 
mensions j gg and j gq . In the large N limit, these are suppressed by a power of 4 in comparison 
to 7 99 and 7 gg , or, equivalently, the corresponding splitting functions are suppressed by a factor 
of 1 — x in the large x limit. Hence, in the large N limit each parton subprocess can be treated 
independently, specifically, each C a b is separately renormalizat ion-group invariant. Because we 
are interested in the behaviour of C a b{N, Q 2 /{x 2 , a s (n 2 )) in the limit N — > oo, we shall henceforth 
omit the parton indices a, b and assume that each subprocess is being treated independently. 
Furthermore, on top of eqs. ()2.3l 12. 7|) the physical process eq. ()2.1|) receives another factorized 
contribution, in which the final-state photon is produced by fragmentation of a primary parton 
produced in the partonic sub-process. However, the cross section for this process is also sup- 
pressed by a factor of 4 in the large N limit. This is due to the fact that the fragmentation 
function carries this suppression, for the same reason why the anomalous dimensions 7 99 and 7 OT 
are suppressed. Therefore, we will disregard the fragmentation contribution. 

Because resummation takes the form of an exponentiation, it is useful to introduce the log 
derivative of the cross section a, i.e., the so-called physical anomalous dimension defined as 



Q 



2 da(N,Q 2 ) 



7 (iV,a s (Q 2 ))cr(iV,g 2 ). 



(2.9) 



The physical anomalous dimension 7 eq. (J2.9|) is independent of factorization scale, and it is 
related to the standard anomalous dimension 7 AP , defined by 



,dF(N,fj, 2 



dfi 2 



7 AP (iV,« s (/. 2 ))F(iV,^), 



(2.10) 



according to 



7 (iV,a s (Q 2 )) 



d\nC(N,Q 2 /fi 2 ,a s (fi 2 )) 



(2.11) 



where we have schematically denoted by F(iV, /z 2 ) the product F^ 1 (N J /j, 2 ) F b H2 (N,/j, 2 ) that ap- 
pears in the cross section for the given parton subprocess, eq. (|2.7|) . In terms of the physical 
anomalous dimension, the cross section can be written as 



a(N,Q 2 ) = K(N;Q 2 } Q 2 )a(N } Q 2 ) = exp E(N;Q 2 } Q 2 ) a(N,Q 2 



(2.12) 



where 



E(N;Q 2 ,Q 2 " 



l{N,a s {k 2 )) 



f k 2 

Q 2 dk 2 
Ql k 2 



(2.13) 



In the large-x limit, the order-n coefficient of the perturbative expansion is dominated by 



terms proportional to 



ln fc (l-a-) 
1-x 



, with k < 2n — 1, that must be resummed to all orders. Upon 
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Mellin transformation, these lead to contributions proportional to powers of In j^. In the sequel, 
we will consider the resummation of these contributions to all logarithmic orders, and disregard 
all contributions to the cross section which are suppressed by powers of 1 — x, i.e., upon Mellin 
transformation, by powers of jr. 

The resummation is performed in two steps, in analogy to ref. |7j. First, we show that the 
origin of the large logs is essentially kinematical: we identify the configurations which contribute 
in the soft limit, we show by explicit computation that large Sudakov logs are produced by the 
phase-space for real emission with the required kinematics as logs of two dimensionful variables, 
and we show that this conclusion is unaffected by virtual corrections. Second, we resum the logs 
of these variables using the renormalization group. 

The Z-th order correction to the leading 0(a s ) partonic process receives contribution from 
the emission of up to / + 1 massless partons with momenta k\, . . . , fcz+i, and 



Pi +P2 = p 7 + fa + . . . h +1 . 

In the partonic center-of-mass frame, where 



Pl = ^(1,0,0,1) 



p a = ^(1,0, 0,-1) 



p 7 = (p ± cosh^ 7 ,pj_,p_L sinhr? 7 ), 



we have 



Q 2 

(Pi +P2 -P-y) 2 = — (1 - Vxcoshr]^) 

X 



l+l 



£ fc°fc°(l-cos%)>0, 



(2.14) 

(2.15) 
(2.16) 

(2.17) 



where 6L- is the angle between ki and kj. Hence, 



x 



Therefore, 



l+i 



1 < cosh?7 7 < —= 



Q 2 

2 



]T A;^ (l-cos%) = ^ r (l-a;) 



O 



:i-x) 



(2.18) 



(2.19) 



Equation ()2.19|) implies that in the soft limit the sum of scalar products of momenta ki of 
emitted partons eq. ()2.17j) must vanish. However, contrary to the case of deep-inelastic scattering 
or Drell-Yan, not all momenta ki of the emitted partons can be soft as x — > 1, because the 3- 
momentum of the photon must be balanced. Assume thus that momenta ki, i = 1, . . . , n; n < l+l 
are soft in the x — > 1 limit, while momenta ki,i> n are non-soft. For the sake of simplicity, we 
relabel non-soft momenta as 



y — h 



1 < j < m + 1; m = I — n. 
The generic kinematic configuration in the x = 1 limit is then 



(2.20) 



fc = 



% = 0; 



1 < i < n 

m+l 

i=i 



P± 



1 < i,j < m + 1. 



(2.21) 



4 



for all n between 1 and /, namely, the configuration where at least one momentum is not soft, 
and the remaining momenta are either collinear to it, or soft. 

With this labeling of the momenta, the phase space can be written as (see the Appendix of 
ref. 0) 

(Pl + Vi\V~ii ki, ■ ■ ■ , k n , hi, ... , k m+1 ) (2.22) 

f s dq 2 fi 2 dk' 2 / / 

= / —d(f) n+1 (p 1 +p 2 ;q,k 1 ,...,k n ) ——d(f) m+1 (k';k' 1 ,...,k' m+1 )d(f)2(q;p^,k'). 

JO 2,71 JO Z7T 

We shall now compute the phase space in the x — > 1 limit. Consider first the two-body phase 
space d<p2 in eq. (|2.22jl . In the rest frame of q we have 



(4tt) £ P 



l-2e 



87iT(l-e) ^ 



sin- 2e ^ 7 d|p 7 | dcosfl 7 5(|p 7 | -P), (2.23) 



where d = 4 — 2e and 



Hence, 



P=f(l-£). (-4) 

Because momenta z < n are soft, up to terms suppressed by powers of 1 — x, the rest frame 
of q is the same as the center-of-mass frame of the incoming partons, in which 

\Pj\ = p± cosily (2.25) 
cos# 7 = tanh?y 7 . (2.26) 

(AttY (Q 2 /4)" e / 2P \ 

dfc(q;p 1 ,k l ) = — — r 7 ^dp ± dr] 1 5 cosh^ 7 = . (2.27) 

87rr(l-e) y/q 2 \ y/Q 2 ) 

The conditions 

IP 

coshr/ 7 = — = > 1; k' 2 > (2.28) 

restrict the integration range to 

Q 2 < q 2 <s (2.29) 
0< k' 2 <g 2 -y / QV- (2.30) 

It is now convenient to define new variables u, v 

q 2 = Q 2 + u(s - Q 2 ) = Q 2 [1 + m(1-x)] + 0((1-x) 2 ) (2.31) 
k' 2 = v{q 2 - y/gV) = Q 2 \uv{\ -x) + 0((1 -x) 2 ) (2.32) 
< m < 1 ; < v < 1, (2.33) 
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in terms of which 



VQ 2 



1 + -u(l - v)(l - x) 







{l-xf 



Thus, the two-body phase space eq. (|2.27|) up to subleading terms is given by 
d(j)2(q;pj,k') 

where 



(Airy (Q 2 /4)- ^ 5( v , - 77+) + jfo - y 



77-t- = In 



87rr(l-e) ^ 



^(l-vXl-x) 



2P 



± 



/ 4 p2 



1 = ±y/u(l-v){l-x). 



[2.34) 



[2.35) 



[2.36) 



KVQ 7 V Q 2 

We now note that the phase-space element d(j) n+ i(pi + p 2 \ q,ki, . . . , k n ) contains in the fi- 
nal state a system with large invariant mass q 2 > Q 2 , plus a collection of n soft partons: 
this same configuration is encountered in the case of Drell-Yan pair production in the limit 
xdy = q 2 /s — > 1, discussed in ref. [IJ. Likewise, the phase space for the set of collinear partons 
d(f) m+ i(k'; k[, . . . , k' m+1 ) is the same as the phase space for deep-inelastic scattering, where the 
invariant mass of the initial state k' 2 vanishes as 1 — x (see eq. ()2.32j) ). We may therefore use 
the results obtained in ref. [Zj: 1 



d(j) n+1 (pi + p 2 ; q,kt,..., kn) = 2tt 
d(j) m+ i(k';k[,...,k' m+1 ) = 2tc 
where N(e) = l/(2(4n) 2 - 2e ) and 



2tt 
2tt 



(q 



s-q 2 ) 



2\2n~l~2ne 



(k 



/2\m— 1- 



me dn '(m) { 



dQ (n \e) (2.37) 



(2.38) 



drt n \e) = dtt 1 ...dtt n C dz n z^- 2)Hn - 1){l - 2t \l -z n f- 2t ... f 1 dz 2 z 1 2 - 2e (l-z 2 ) 1 - 2 2.39) 

Jo Jo 



dfi'M(e) = dn[ ...dn' m [ dz'J^-^{\ - z' m r 



-2c 



dz' 2 z 2 e (l - : 



I \l-2e 
2 



(2.40) 



The definition of the variables z[ is irrelevant here and can be found in ref. [7]. 

Equations ()2.31I2.32|) imply that the phase space depends on (1 — x)~ e through the two 
variables 



[s - <r 



k' 2 ocg 2 (i-x) 

" 2 ocQ 2 (l-a;) 2 , 



(2.41) 
(2.42) 



where the coefficients of proportionality are dimensionless and x-independent. By explicitly 
combining the two-body phase space eq. (|2.35|) and the phase spaces for soft radiation eq. (|2.37|) 

^dn the case of deep-inelastic scattering, in ref. [Jj one of the outgoing parton momenta {k' m+1 , say) was identi- 
fied with the momentum of the outgoing quark and called p', hence eq. (|2.38|) is obtained from the corresponding 
result in ref. [xj by the replacement p' — > k' m+1 . 
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and for collinear radiation eq. (|2.38j) we get 



(-[ T \2n+m— (2n+m)e 

^ m+2 (p 1 +p 2] p 1 ^^...^ n ^[ ) ...^U 1 ) = m n+m - {n+m+l)e -- 



-m+me 



(16tt) 



-l+g r at/ A 1 n+m 



r i 



N(e) 
2tt 



dr]^dn {n \e) dn' {m) ( 



du K —^J / dv [5( Vl - V+ ) + 8{ Vl - V -)]. (2.43) 

u Jo vl — v 



In the limiting cases n = and m = we have 

dMPi+P2;q) = 2ir5(s-q 2 ) = Q2{1 _ %) 6(1 -u) (2.44) 

dMk';p') = 2n5(k' 2 ) = - <fo); (2.45) 

the corresponding expressions for the phase space are therefore obtained by simply replacing 

(1 - u)- 1 dfl^Xe) ^ 8(1 - u); v~ x d£l^ m \e) -> <S(v) (2.46) 

in eq. (|2.43j) for n = 0, m = respectively. 

The logarithmic dependence of the four-dimensional cross section on 1 — x is due to inter- 
ference between powers of (1 — x)~ e and - poles in the d-dimensional cross section. Hence, we 
must classify the dependence of the cross section on powers of (1 — x)~ e . We have established 
that in the phase space each real emission produces a factor of [Q 2 (l — x) 2 ] 1 e if the emission is 
soft and a factor of [Q 2 (l — x)] 1 6 if the emission is collinear. The squared amplitude can only 
depend on (1 — x)~ e because of loop integrations. This dependence for a generic proper Feynman 
diagram G will in general appear [2j through a coefficient 

A G (P E ) = [D G (P E )} dL/2 - 1 , (2.47) 

where L and / are respectively the number of loops and internal lines in G, and Dg{Pe) is a 
linear combination of all scalar products Pe of external momenta. In the soft limit it is easy to 
see, by manipulations analogous to eq. (j2.17j) . that all scalar products which vanish as x — > 1 
are either proportional to Q 2 (l — x) or to Q 2 (l — x) 2 . Equation ()2.47|) then implies that each 
loop integration can carry at most a factor of [Q 2 (l — x) 2 ] 6 or [Q 2 (l — x)] e . 

This then proves that the perturbative expansion of the bare coefficient function takes the 
form 

oo 

C (0) (x,Q 2 ,« ,e) = aa (Q 2 y e Y, a oCl°\x,Q 2 ,e) (2.48) 

1=0 

C} 0) (x,Q 2 ,a Q ,e) = —^p^== £ £ c£,(e)(l - x)- 2ke ~ k \ (2.49) 
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where the coefficients Cj^., have poles in e = up to order 21, and the factor l/r(l/2) was 
introduced for later convenience. Terms with k + k' < I at order a! s are present in general 
because of loops. The Mellin transform of eq. (|2.48j) can be performed using 



dxx 



N-li 



x)' 



-l/2-2fce-fc'e 



with the result 

C^(N,Q 2 ,a ,e 
aa (Q 



2^— e oo I l—k 

-EEE4 ^) 

Z=0 k=0 k'=0 







k 




_{%) 


a 







y 



(Q 



2\-e. 



a 



l-k-k' 



+ 



[2.50) 

(2.51) 
1 

N 



3 Resummation from renormalization group improvement 

Equation (|2.51|l shows that indeed as iV — - > oo, up to corrections, the coefficient function 

depends on N through the two dimensionful variables j& and ^r. The argument henceforth is 
a rerun of that of ref. [7j, in this somewhat more general situation. The argument is based on 
the observation that, because of collinear factorization, the physical anomalous dimension 



j(N, a s (Q 2 )) = Q 2 ^ \nC(N, Q 2 //i 2 , a s ^ 2 )) 



(3.1) 



is renormalization-group invariant and finite when expressed in terms of the renormalized cou- 
pling a s (/i 2 ), related to «o by 



« (/i 2 , «,(/i 2 )) = /i 2e « s (/i 2 ) Z^\a s {^ 2 ), e), 



(3.2) 



where Z^ as \a s (/j, 2 ), e) is a power series in a s (/i 2 ). Because «o is manifestly independent of /i 2 , 
eq. (|3.2|) implies that the dimensionless combination (Q 2 )~ e a (a s (ii 2 ), /i 2 ) can depend on Q 2 only 
through a s (Q 2 ): 

(Q 2 r*a (a s ( f i 2 ), f i 2 ) = a s (Q 2 )Z^\a s (Q 2 ),e). (3.3) 

Using eq. (|3.3j) in eq. (|2.51|) . the coefficient function and consequently the physical anomalous 
dimension are seen to be given by a power series in a s {Q 2 ), a s (Q 2 /N) and a s (Q 2 /N 2 ): 



oo oo oo 



7(iV,a s (Q 2 ), e ) = E E E 7™„ P (e) <(Q 2 ) <(Q 2 /N 2 ) a p s {Q 2 /N). 



(3.4) 



m=0 n=0 p=0 



Even though the anomalous dimension is finite as e — > 0, the individual terms in the expansion 
eq. ()3.4|) are not separately finite. However, if we separate iV-dependent and iV-independent 
terms in eq. (|3.4j) : 



j(N, a s (Q 2 ),e) = 7 (c) (« s (Q 2 ), e) + 7 W (iV, a s (Q 2 ), e ) 



(0, 



(3.5) 
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we note that the two functions 

^(a s (Q 2 ),e) ee 7 (c) («s(Q 2 ),e)+7 (0 (l,« s (Q 2 ),6) (3.6) 
^(N,a s (Q 2 ),e) = ^(N,a s (Q 2 ),e)-^\l,a s (Q 2 ),e) (3.7) 

must be separately finite, because 

7 (iV, a s (Q 2 ), e) = 7 (c) («s(Q 2 ), 6) + 7 (0 (iV, a s (Q 2 ), e), (3.8) 

is finite for all N, and 7") vanishes for N = 1. 
We can rewrite conveniently 

7 (/) (^,« s (Q 2 ),e)= / — g(a s (Q 2 ),a s (Q 2 /n 2 ),a s (Q 2 /n),e), (3.9) 
Ji n 

where 

s(a s (Q 2 ), a,(g 2 /ri 2 ), a s (Q 2 /n), e) = n^ l \n, a s (Q 2 ), e). (3.10) 

is a Taylor series in its arguments whose coefficients remain finite as e — > 0. In four dimension 
we have thus 

7 (iWQ 2 )) = ^ l) (N,a s (Q 2 ),0)+^\N,a s (Q 2 ),0) + O^ 
= ^ l \N,a s (Q 2 ),0) + O(N°) 

= f N -g(a s (Q 2 ),a s (Q 2 /n 2 ),a s (Q 2 /n)) + 0(N°), (3.11) 

where g(a s (Q 2 ),a s (Q 2 /n 2 ),a s (Q 2 /n)) = lim e ^ g(a s (Q 2 ),a s (Q 2 /n 2 ),a s (Q 2 /n),e) is a generic 
Taylor series of its arguments. 

Renormalization group invariance thus implies that the physical anomalous dimension 7 
eq. ()3.1|) depends on its three arguments Q 2 , Q 2 /N and Q 2 /N 2 only through a s . Clearly, any 
function of Q 2 and iV can be expressed as a function of a s (Q 2 ) and a s (Q 2 /N) or a s (Q 2 /N 2 ). 
The nontrivial statement, which endows eq. ()3.11|) with predictive power, is that the log deriva- 
tive of 7, g(oi s (Q 2 ),a s {Q 2 /n 2 ),a s (Q 2 /n)) eq. (|3.1U|) . is analytic in its three arguments. This 
immediately implies that when 7 is computed at (fixed) order a k s , it is a polynomial in In of 
fc-th order at most. 

In order to discuss the factorization properties of our result we write the function g as 

g(a s (Q 2 ) } a s (Q 2 /n 2 ) } a s (Q 2 /n)) = gi (a s (Q 2 ), a s (Q 2 /n)) + g 2 (a s (Q 2 ) } a s (Q 2 /n 2 )) 

+g s (a s (Q 2 ),a s (Q 2 /n),a s (Q 2 /n 2 )) 

00 00 

gMs(Q 2 )^ s (Q 2 /n)) = Y.Y.9mv P <(Q 2 )a p AQ 2 h) (3.12) 

00 00 

g 2 (a s (Q 2 ),a s (Q 2 /n 2 )) = ]T £ g mn0 a^(Q 2 ) a n s (Q 2 /n 2 ) 



m=0 n=l 

00 00 00 



g 3 (a s (Q 2 ),a s (Q 2 /n),a s (Q 2 /n 2 )) = Y,Y,Y,9mn P <{Q 2 )a n s {Q 2 /n 2 )al{Q 2 /n). 

m=0 n=l p=l 



The dependence on the resummation variables is fully factorized if the bare coefficient functions 
has the factorized structure 

C(°\N, Q\ a , e) = C^ C \Q\ a , e) C[°' 1 \q 2 /N, a , e) C^' l \Q 2 /N 2 , a , e). (3.13) 

This is argued to be the case in the approach of refs. [HI Ell- If this happens, the resummed 
anomalous dimension is given by eq. (|3.11|) with all g mnp = except go n o,goo P '- 

1 (N,a s (Q 2 )) = — gi ^a s (Q^n))+ — g 2 (0, a s (Q 2 /n 2 )). (3.14) 
Ji n Ji n 

Note that because the coefficient function depends on the parton subprocess (compare eq. (|2.3Jl ) 
the factorization eq. (|3.13|) applies to the coefficient function corresponding to each subprocess, 
and the decomposition eq. ()3.14j) to the physical anomalous dimension computed from each of 
these coefficient functions. 

A weaker form of factorization is obtained assuming that in the soft limit the iV-dependent 
and iV-independent parts of the coefficient function factorize: 

C (0 \N, Q 2 , e) = C^ C \Q 2 , oo, e) C (0 > l \Q 2 /N 2 , Q 2 /N, a , e). (3.15) 

This condition turns out to be satisfied [7J in Drell-Yan and deep-inelastic scattering to order a 2 . 
It holds in QED to all orders ^3] as a consequence of the fact that each emission in the soft limit 
can be described by a universal (eikonal) factors, independent of the underlying diagram. This 
eikonal structure of Sudakov radiation has been argued in refs. jSHHI to apply also to QCD. If 
the factorized form eq. (|3.15J) holds, the coefficients g mnp eq. ()3.12j) vanish for all m^O, and the 
physical anomalous dimension takes the form 

7 (iV,a s (Q 2 )) = / — gi ^a s (Q 2 /n))+ / - g 2 (0, a s (Q 2 /n 2 )) 
Ji n Ji n 

+ / —g 3 (0,a s (Q 2 /n 2 ),a s (Q 2 /n)). (3.16) 
Ji n 

It is interesting to observe that in the approach of refs. [HIE] for processes where more than 
one colour structure contributes to the cross-section, the factorization eq. ()3.13|) of the coefficient 
function is argued to take place separately for each colour structure. This means that in such 
case the exponentiation takes place for each colour structure independently, i.e. the resummed 
cross section for each parton subprocess is in turn expressed as a sum of factorized terms of the 
form of eq. ()3.13|) . This happens for instance in the case of heavy quark production ^OJ IT2"] . 
In prompt photon production different colour structures appear for the gluon-gluon subprocess 
which starts at next-to-next-to-leading order, hence their separated exponentiation would be 
relevant for next-to-next-to-leading log resummed results. 

When several colour structures contribute to a given parton subprocess, the coefficients of 
the perturbative expansion eq. (|2.51|) for that process take the form 

CtW) = C\ll}{e)+C\ll S {e\ (3.17) 
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(assuming for definiteriess that a colour singlet and octet contribution are present) so that the 
coefficient function can be written as a sum = C{ + Cg \ The argument which leads from 
eq. (j2.51|) to the resummed result eq. ()3.11|) then implies that exponentiation takes place for each 
colour structure independently if and only if 7j = d In /d In Q 2 and 7s = d In jd In Q 2 
are separately finite. 

This, however, is clearly a more restrictive assumption than that under which we have derived 
the result eq. (jH.lljl . namely that the full anomalous dimension 7 is finite. It follows that 
exponentiation of each colour structure must be a special case of our result. However, this 
can only be true if the coefficients g^ k of the expansion eq. (j3.12|) of the physical anomalous 
dimension satisfy suitable relations. In particular, at the leading log level, it is easy to see that 
exponentiation of each colour structure is compatible with exponentiation of their sum only 
if the leading order coefficients are the same for the given colour structures: g^ m = g^ 01 and 
#010 = #010- This is indeed the case for heavy quark production (where (7001 = 0). 

Note that, however, if the factorization holds for each colour structure separately it will not 
apply to the sum of colour structures. For instance, the weaker form of factorization eq. (j3.15|) 
requires that C^,(e) = F fc+fc /(e)G;_ fe _ fc /(e), but 

Fl +k ,{z)GU-y{e) + F* +kl (e)Gl k _ h ,(e) + F k+k ,(z)Gi- k - k >(e). (3.18) 

Hence, our result eq. ()3.11|) for the sum of colour structures is more general than the sepa- 
rate exponentiation of individual colour structures, but it leads to results which have weaker 
factorization properties. 

The resummation coefficients g m np can be determined by comparing the expansion of the 
resummed anomalous dimension 7 in powers of a s (Q 2 ) with a fixed-order calculation: 

^min i 1 

lFO (N,a s ) = £a£ £7i W +0 (a k r^ 1 ) + 0(N°), (3.19) 

i=l j=l 

where 7fo(^, «s) is the physical anomalous dimension for the same individual partonic subpro- 
cess (recall eq. (j2.3|0 . Clearly, if the more restrictive factorized forms eq. ()3.13|) or eq. ()3.15|) 
hold, a smaller number of coefficients determines the resummed result, and thus a lower fixed- 
order calculation is sufficient to predict higher-order logarithmic terms than if the more general 
eq. ()3.12|) is used. Conversely, a higher fixed-order calculation can be used to verify if the strong 
factorization eq. ()3.13j) holds as advocated in refs. (HI EH], or whether one must use the less 
predictive but more general result eq. ()3.12j) that we have derived. 

Once the resummed physical anomalous dimension has been determined, the resummed cross 
section can be obtained from it using eq. ()2.12|) . with a factorization scheme choice which specifies 
the way it is split into its two components eq. (j2.11|) . Commonly used choices are the physical 
scheme choice, in which C = 1 so 7 = 7 AP , or the MS scheme, in which the unresummed and 
resummed forms of the anomalous dimensions 7 AP coincide. An explicit construction of the 
relation between physical anomalous dimension and resummed cross section, and the matching 
between resummed and unresummed results, can be found in Section 6 of ref. [7j. 
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4 The structure of resummed results 



We determine the predictive power of each resummed result by means of the following strategy. 
First, we assume that the coefficients g m np needed for N fc_2 LL resummation have already been 
determined. Next, we identify the coefficients that are needed to extend the accuracy to N fc_1 LL, 
and we write a system of equations that fix them in terms of the known coefficients, and of the 
Yj of the fixed-order expansion. The rank of this system of equations determines the minimum 
order k min in a s of a fixed order computation which is needed to fix the N fc_1 LL resummation. 
This means that at any higher fixed order / > k min} the coefficients of all powers of (in ~J with 
k < n < f are then predicted by the resummed formula. 

The general structure of the anomalous dimension resummed to N fc_1 LL accuracy is 

1 (N,a s (Q 2 )) = J2l P (N,a s (Q 2 )) (4.1) 
P =i 

lp (N } a s (Q 2 )) = J2T,9^-j^(Q 2 ) / -ai(Q 2 /n 2 )ar- J (Q 2 /n). (4.2) 

i=0,7=0 Jl 71 

At the N fc_1 LL order, in each term •y p (N,a s (Q 2 )) the coupling constant a s (Q 2 /n a ) can be ex- 
panded in powers of a s (Q 2 ) using the N fc ~ p LL solution of the renormalization group equation 

/i 2 ^ = -/?o« s 2 -/?i« s 3 + ... , (4.3) 

because subsequent orders would lead to N fc LL contributions to 7. In particular, the leading log 
expression 

*s{Q 2 /n a ) = — "fffi 1 (4.4) 
1 + aa s {Q 2 )(3 m- 

is sufficient for the computation of 7fc(iV, a s (Q 2 )). With a s (Q 2 /n a ) given by eq. (|4.4jl one gets 
/ - a\{Q 2 /n 2 ) ai(Q 2 /n) = £ C$fi ffl a s (Q 2 ) i+ ^ m ln m+1 - (4.5) 

Jl U m=0 iV 



(note that y U ^ J — 1 for n — — 1 and otherwise), and therefore 

k— 1 k—i 00 1 

lk (N, a s (Q 2 )) = ££ QHk-i-i E C^- l - 3) Po*s{Q 2 ) k+m ln m+1 W (4.7) 

i=0 j=0 m=0 iV 

Let us consider first the leading log case, k — 1. In this case, g is linear in a s and therefore 
eq. (|H.14|) always holds, i.e. there is no distinction between factorized and unfactorized results. 
The anomalous dimension has the form 

7i(iV,a s (Q 2 )) = £ (#001 +2 m g 0W ) /W+W) ln m+1 - + NLL. (4.8) 

m=0 171 ' 1 iV 
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Comparing with the fixed order expansion, we find 

m = : 7} = -(#001 + 9ow) (4.9) 

m=l: 72 = y (<?ooi + 2<?oio)- (4.10) 

These two conditions determine g oi an d #010 : leading-log resummation requires computing 7 
to order a 2 . Note that for DIS and Drell-Yan an 0(a s ) computation is instead sufficient [Zj, 
because only one scale is present and thus only one coefficient has to be determined. 

Let us now examine higher logarithmic orders, by discussing the various factorizations in 
turn. Assume first the validity of the most restrictive result eq. (|3.14|) . where the dependence on 
the three scales Q 2 , Q 2 /N and Q 2 /N 2 is fully factorized. In such case, the anomalous dimension 
is just the sum of a function of a s (Q 2 /N) and a function of a s (Q 2 /N 2 ), and all coefficients 
Qmnp vanish except g 00p , g 0p Q. In this case, A^ = 2k coefficients are required for the N fc_1 LL 
resummation. According to the strategy outlined above, we now assume that the coefficients 
goop, gopo with p < k — 1, relevant for N fc_2 LL resummation, have already been determined. The 
two extra coefficients goofc, 9ow appear in jk(N, a s ), whose explicit form is given by eq. ()4.7|) with 
only gook,gow different from zero: 

00 1 
lk (N,a s (Q 2 )) = T,(C^ k) gook + C^g oko )^a k s +m (Q 2 )\n m+1 -. (4.11) 

m=0 iV 

Other terms of order a!s +m ln m+1 (l/7V) are generated by expanding the coupling a s (Q 2 /n a ) 
in 7x,...,7fc_i up to N fc_1 LL; however, these are fully determined by the known coefficients 
9ooi,9oio,i < k — 1 (and by the coefficients of the beta function fi , . . . ,(3k-i). Equating terms 
of order a k+m \n m+1 [1 / N) in the fixed-order and resummed expressions of 7 we get the set of 
equations 

m = : 7l fc = -(#00* + 9oko) + Fl k \g m , g m ; l<i<k-l) (4.12) 

m = 1 : 7 2 +1 = (2oofc + 2g oko ) ^ + F^ k \g 00l , g m - l<i<k — l) (4.13) 
. . . , 

where Fj^ are known functions of the coefficients g which we have assumed to be already known. 

Hence, the two terms m — 0, 1 provide two independent conditions that fix g ofc and gofco in 
terms of #00;, goio, i < k — 1. The same procedure can be repeated for p = k — 1, k — 2, . . . , 1; 
at each step, goo p and gopo are computed as functions of 7i,72 +1 and gooi, gmo, < i < p — 1. 
We conclude that in the case of eq. ()3.14|) the coefficients gooi, goio] < i < k, relevant for 
N fc_1 LL resummation, are obtained from the fixed-order expansion of 7 up to order k m i n = k + 1 
(corresponding to m = 1 in eq. (|4.11|) ). This means that even though at each extra logarithmic 
order two new coefficients appear, a single extra fixed order in a s is sufficient to determine both 
of them. 

Consider for example the next-to-leading log resummation. In our approach, this requires 
the calculation of the physical anomalous dimension up to order a^. Explicitly, 

7i(Af,a s (Q 2 )) = -(#001 + 9oio) <y s {Q 2 ) In — 
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+ ^(^001 + 2g io) a 2 s (Q 2 ) In 2 -J- - ^(#001 + 9ow) a 2 (Q 2 ) In 1 



N Po 



N 



01 



3ft 



~{9oai + 4<7„io) a'(Q 2 ) In 3 ^ + -y^ooi + 2<foio) « S 3 (Q 2 ) In 2 ^4.14) 



72 (iV,«,(Q 2 )) 



-(#002 + 6-020) a 2 s {Q 2 ) In — + /3 q (5'qo2 + 2sf 020 ) a^(Q 2 ) In 2 — . 



(4.15) 



The leading log coefficients are fixed by the linear and quadratic In terms in 71, eqs. (|4.9|4.1(J|) . 
and the NLL coefficients by the the linear and quadratic In ■h terms in 72: 



#002 + ^020 - 

6002 + 2(^020 



~ 7T (fi'ooi + fltoio) 



00 



72 



3ft 
2ft 



J001 



2g ( 



010 > 



(4.16) 
(4.17) 



All other logarithmically enhanced contributions to the order anomalous dimension are pre- 
dicted, except the In — which depends on the N 2 LL coefficients. The same pattern continues 
at higher orders. 

It is interesting to observe that even assuming the fully factorized form of the coefficient 
function eq. (|3.13j) . our results are still less restrictive and thus less predictive than those of 
refs. [HI E3, where one less fixed order is required to determine the resummed result. For 
instance, at NLL, using the notation of ref. ^U|j the resummation coefficients are given by 



#010 
#001 
6020 
#002 



AO) + A® 



A 



(i) 



7T 



.4 



(i) 



r4 2) + 4 2) 



(4.18) 
(4.19) 



.4 



(2) 



ft(4 x) + 4 1} - 4 1J ) 



l(l)> 



7T 



7T 



m2-2 TE ft^ 



A 



(i) 



7T 



^(4.20) 



.4 



(2) 



7T^ 



a R 1 4(1) 

- IePo- 



2tt 



7T 



(4.21) 



where A$ is the coefficient of ln(l/AT) in the Mellin transform of the P aa Altarelli-Parisi splitting 
function at order a\, 7^ is the Euler constant, and B?' is a constant to be determined from the 



comparison with the fixed-order calculation. In eq. 



a, b are the incoming partons (on 



which C implicitly depends), and d is the outgoing parton in the leading order process when 
incoming partons a and b (which is uniquely determined by a and b). 

Thus, in this approach go«o is entirely determined in terms of the 0(a l s ) coefficient of the 
ln(l/AQ term in the anomalous dimension, and only goo« must be determined by comparison to 
the fixed-order calculation: at the LL level, the resummed result is obtained from the knowledge 
of 7fo to order a s , at the NLL level to order a 2 and so on. This means that to LL, the coefficient 
7! eq. (|4.1()jl is in fact predicted by eqs. (j4.18J4.19j) in terms of the coefficients of the Altarelli- 
Parisi splitting functions. This prediction is borne out by the explicit 0(a 2 ) calculation of the 
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prompt photon production cross section ^3]. At the NLL level, the coefficient 7 f of ln 2 (l/iV), 
is predicted: 



-Po 



2(4 2) + 4 2) )-A 



I (2) 



7T 



/3 (21n2 + 4 7E )^ 



^ + A 



(i) 



7T 



(4.22) 



+/3 (21n2 + 3 7B ) 



.4 



(i) 



A>3 



d 



7T 



2tt 



M 2 



2(4 1) + 4 1) )-4 ] 



7T 



The correctness of this result could be tested by an order calculation. If it were to fail, the 
more general resummation formula with (7020 determined by eq. (|4.17|) should be used, or one of 
the resummations which do not assume the factorization eq. (|3.13|) . 

Assume now that the weaker factorization eq. (|3.15j) applies. In this case, only g 0ng ; 1 < 
n + q < k are nonzero. This amounts to keeping only the term i = in the general expression 
eq. (J4~2"j) for 7 P (iV, a s (Q 2 )). The total number of coefficients at N fc_1 LL is 



p=i z 



(4.23) 



In order to improve the accuracy from N fe ~ 2 LL to N fc_1 LL, k + 1 new coefficients are needed, 
namely goik-i for i — 0, . . . , k. As in the previous case, we isolate the N fc_1 LL terms that depend 
on the new coefficients. All such terms are contained in 7^, which in this case is given by the 
term % = in eq. ()4.7|) : 



j=0 m=0 iV 

The first k + 1 terms in the sum over m provide the set of conditions 



(4.24) 



m = : 



m = 1 : 



m = k : 



where again are known functions of the coefficients g which we have assumed to be already 
known. Equations (j4.25H4.27j) are linearly independent (see the Appendix for an explicit proof), 
and therefore determine goik-f,0 < i < k in terms of gojk-j',0 < j < k — 1. It follows that 
a computation of 7 up to order fc min = 2k is necessary for the determination of the k ( k ^ 
coefficients needed for N fc_1 LL resummation: even though the number of coefficients which 
must be determined grows quadratically according to eq. ()4.23|) . the required order in a s of the 
computation which determines them grows only linearly. 



li = E /AIM- + g? W-i; < 3 < k - 1) 

j=0 



7 2 fc+1 



E gojk-jC?' 3) fa + G [ 2 ] {gojk~j] < j < k- 1) 
3=0 



iT + i = E 9o jk -fir~ 0) P k + GUxigojk-j; < j < k - 1), 
3=0 



(4.25) 
(4.26) 

(4.27) 



15 



Note that in fact a similar argument can be applied to the cases of deep-inelastic scattering 
and Drell-Yan in the absence of factorization, discussed in ref. [7j. In that reference it was 
observed that k(k + l)/2 coefficients must be known in order to determine the N fc_1 LL result, 
and thus a fixed N fc ^ +1 - )/ ' 2_1 LO computation is sufficient to determine them. However, using 
the argument given here it is easy to see that in fact a calculation of 7 up to fixed order 
&mm = 2k — 1 only is sufficient to determine all these coefficients. Consider for example the case 
of deep-inelastic scattering. To N fc_1 LL one finds 

7 (iV, a s (Q 2 )) = £ lp (N, a s (Q 2 )) (4.28) 
p=i 

p " 1 dn 



7p (iV, as (g 2 )) = ^^_ ia ' s (Q 2 ) / -ar(Q 2 /n). (4.29) 



In order to determine the k coefficients gok, ■ ■ ■ ,9k-u that are needed to improve the accuracy 
from N fc ~ 2 LL to N fc_1 LL, we proceed as in the previous case: the new coefficients appear only in 

fc— 1 00 1 

7fc (iV, a s (Q 2 )) = £ 9ik-i E C ^~ l) ffl «s(Q 2 ) fe+m ln m+1 -. (4.30) 

i=0 m=0 iV 

Each term in the sum over m in eq. (|4.3U|) provides an independent condition on the coefficients 
gij (the linear independence of these conditions is straightforwardly proved in the Appendix). 
Hence, in order to determine all Qik-i] < i < k — litis sufficient to determine all terms up to 
m = k — 1 in 7^, i.e. compute 7 up to order a 2k ~ l . The same happens in the case of Drell-Yan, 
which is obtained by replacing C^' p ~^ with C^~ % ^ in eq. (|4.30p . 

Let us now consider the most general case, in which the coefficient function does not satisfy 
any factorization property. Then, to N fc_1 LL, the anomalous dimension eqs. (J4.1I4.2)) depends 
on 

Nt = j 2 P_ip ± S 1= k(k + l)(k + 5) (431) 
P =i ^ 6 

coefficients overall; out of these, the new ones which must be determined in order to go from 
N fc - 2 LL to N fc_1 LL are the k(k + 3)/2 coefficients 

gijk-i-j] i = 0, . . . , k - 1; j = 0, . . . , k - i. (4.32) 

These new coefficients are contained in 7^, now given by its general expression eq. (|4.7jl . and each 
term with fixed m in the expansion of 7^ eq. ()4.7|) provides a new condition on these coefficient. 
However, these conditions are not linearly independent for any choice of m: rather, the rank of 
the matrix which gives the linear combination of coefficients eq. (J4.32)) to be determined turns 
out to be 2k (see the Appendix). This means that the N fc_1 LL order resummed result depends 
only on 2k independent linear combinations of the k(k + 3)/2 coefficients eq. (|4.32j) . Because 
a term with fixed m in 7^ is of order a k s +m , this implies that a computation of the anomalous 
dimension up to fixed order /c min = 3k — 1 is sufficient for the N fc_1 LL resummation. Note that 
when going from N fe_1 LL to N fc LL and 7^ is now determined at this higher order, in general 
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DIS, DY 




eq. (gUD 


eq. (f3~Tnj) 


eq. dHU) 


refs. [2, 3 


ref. [Z| 




2k 


k(k+3) 
2 


fc(fc+l)(fe+5) 
6 


k 


k(k+i) 

2 


h ■ 


k + 1 


2A; 


3A; - 1 


k 


2A; - 1 



Table V. Number of coefficients and minimum order of the required perturbative calculation 
kmm for different versions of the N fc_1 LL resummation. 



some new linear combinations of the k(k + 3)/2 coefficients eq. ()4.32j) will appear through terms 
depending on fii. Hence, some of the combinations of coefficients that were left undetermined in 
the N fc_1 LL resummation will now become determined. However, this does not affect the value 
k min of the fixed-order accuracy needed to push the resummed accuracy at one extra order. In 
conclusion, even in the absence of any factorization, despite the fact that now the number of 
coefficients which must be determined grows cubically according to eq. (|4.31|) . the required order 
in a s of the computation which determines them grows only linearly. 

The number of coefficients Nj. that must be determined at each logarithmic order, and the 
minimum fixed order which is necessary in order to determine them are summarized in Table 
1, according to whether the coefficient function is fully factorized [eq. (|3.14|) ]. or has factorized 
iV-dependent and iV-independent terms [eq. f!3.16|) ]. or not factorized at all [eq. (|3.11|) ]. In 
the approach of refs. [HI EI| the coefficient function is fully factorized, and furthermore some 
resummation coefficients are related to universal coefficients of Altarelli-Parisi splitting functions, 
so that /c m in = k. For completeness, we also list in the table the results for DIS and Drell-Yan, 
according to whether the coefficient function has factorized iV-dependent and iV-independent 
terms (as in refs. [3 El) or no factorization properties (as in ref. [7]). Current fixed-order results 
support factorization for Drell-Yan and DIS only to the lowest nontrivial order 0(a 2 s ). For 
prompt-photon production, available results do not allow to test factorization, and test relation 
of resummation coefficients to Altarelli-Parisi coefficients only to lowest 0(a s ). 

5 Conclusion 

In this paper, we have presented a generalization to prompt photon production of the approach 
to Sudakov resummation which was introduced in ref. [7j for deep-inelastic scattering and Drell- 
Yan production. The advantage of this approach is that it does not rely on factorization of 
the physical cross section, and in fact it simply follows from general kinematic properties of the 
phase space. It is interesting to see that this remains true even with the more intricate two-scale 
kinematics that characterizes prompt photon production in the soft limit, especially in view of 
the fact that the theoretical status of Sudakov resummation for prompt photon production is 
rather less satisfactory than for DIS or Drell-Yan. Also, this approach does not require a separate 
treatment of individual colour structures when more than one colour structure contributes to 
fixed order results. 

The resummation formulae obtained here turn out to be less predictive than previous re- 
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suits |EJ EI] : a higher fixed-order computation is required in order to determine the resummed 
result. This depends on the fact that here neither specific factorization properties of the cross 
section in the soft limit is assumed, nor that soft emission satisfies eikonal-like relations which 
allow one to determine some of the resummation coefficients in terms of universal properties 
of collinear radiation. Currently, fixed-order results are only available up to 0{af) for prompt 
photon production. An order a z s computation is required to check nontrivial properties of the 
structure of resummation: for example, factorization, whose effects only appear at the next-to- 
leading log level, can only be tested at O(af). The greater flexibility of the approach presented 
here would turn out to be necessary if the prediction obtained using the more restrictive resum- 
mation of refs. jHl IE| were to fail at order af. 



A Appendix 

In this Appendix, we prove some properties of matrices built from the coefficients eq. (|4.6|h 
which appear in the perturbative expansion of the highest order contribution 7*. to the resummed 
anomalous dimension eq. (j4.H4.2j) . 

1. The k x k matrix 

A^ = CS' k ~ i) ] 0<m<k-l, 0<i<k-l (A.l) 

is non-singular. 

Proof: From eq. (j4.fi)) we see that A^j is a degree-m polynomial in i: 

I 1 \m+l m 

(m+1)! ^ 
It follows that a generic linear combination of the rows of A^ 

k—l k—l m k—l I 

X m A mi = X m Ail = A\ % X m (^-3) 

m=0 m=0 1=0 1=0 rn=0 

can only vanish if x m = for all m. 

It follows that each term in the sum over m in eq. ()4.30|) provides a linearly independent 
condition on the coefficients gik-i, < i < k — 1. 

2. The k x k matrix 

A^ = Ci fc ~ i,0) ; 0<m<A;-l, ®<i<k-l (A.4) 

is non-singular. 

Proof: This statement follows immediately from the previous one, because it is easy to 
show that 

Wfc-i,0) _ 9m W0,fc-i) /a p-\ 
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3. The (jfe + 1) X (jfe + 1) matrix 

B%}=C£ k -» ] 0<m<k, 0<j<k (A.6) 

is non-singular. 

Proof: This statement can be proved by induction on k. For k — 1 we have 

= ( T/2 V ) < A ' 7 > 

which is manifestly non-singular. We now assume that B^ -1 ' is non-singular, and we 
consider a linear combination of the columns of 

£ Xj b£> = (xo + 2 m x k ) + J2 x s (A.8) 

3=0 j=l 



where we have used eq. (|A.5|) . For i,j > 1 the following identity holds: 

c m = acgj-i) _ c g-ij) i (A . 9) 

Equation (|A.9J) can be verified directly, using the standard properties of the binomial 
coefficients 



n — 1 \ _ f n \ I n — 1 \ / n 

£;- 1 ] = I £; ] ~ I fc J ' I k 



for n < k. (A. 10) 



Using eq. (jA.9|) we get 



fc fc-l k-l 

= (x + 2 m Xfc) C^'*) + 2^Xj C^ k j 1 - ) — ^ x j &m 1,k ^ 

3=0 j=l J=l 

= (o; Q + 2 m x k ) + 2 J2 Xj - £ x j+1 Cg**" 1 ^ 

i=l i=0 

= {x Q + 2 m x k )C^ + Y j x 3 Bt- 1 \ (A.li; 

i=o 



where 



Xo — — X\ 

Xj = 2xj - x j+1 ; 1 < j < k - 2 (A. 12) 

x fc _i = 2xk-i- 

The linear combination in eq. (jA.ll|) can only vanish if the two terms are separately zero, 
since C^ ,k ^ is a degree-(A; — 1) polynomial in m, while B^ ^ is at most of degree k — 2. 
Hence, for eq. (jA.ll|) to vanish, it must be 

x + 2 m x fc = (A.13) 
ExyBSr^O, (A.14) 

3=0 

19 



and therefore 

x = 0; Xk = 0; Xj = (A. 15) 

by the induction hypothesis. Using eqs. (jA.12|) this gives 

xo — x\ = . . . = Xk = 0. (A. 16) 

It follows that eqs. (j4.25H4.27j) provides a linearly independent condition on the coefficients 
flto-i, < i < k. 

4. Define an M x fc ( fc + 3 ) matrix D^, whose columns are the M-component vectors 

D<J) = Cg** - *"^; < z < A; - 1; 0<j<k-i; < m < M. (A.17) 
The rank (number of linearly- independent columns) of is 2fc. 
Proof: We use induction on fc. For k = 1, is a 2 x 2 matrix with columns 

( _i \m+l 

^ = K '^ ^' 0) ) = a, n, (A.i8) 

that are linearly independent; the rank of is 2. Let us check explicitly also the case 
k = 2. In this case g 

r>(2) _ /wo.1) mi,q) r(°' 2 ) r^ 1 ' 1 ) r( 2 - ^ (a itn 

The first two columns are the same as in the case k = 1: they span a 2-dimensional 
subspace. The last three columns are independent as a consequence of statement El of 
this Appendix. Furthermore, C^' 2 -* and C^ 2 ' -* = 2 m C^' 2 ^ are independent of all other 
columns, because they are the only ones that are proportional to a degree- 1 polynomial 
in m. Finally, C^ 1 ' is a linear combination of the first two columns, as a consequence of 
eq. (TP1) with i = j = 1. Thus, the rank of £> (2) is 2 + 2 = 4. 

We now assume that has rank 2{k — 1), and we write the columns of as 

D£> = (Cg-*- 1 ^, Cg*^) (A.20) 
< i < A: — 2, < j < A; - 1 - z 0<l <k. (A.21) 

By the induction hypothesis, only 2(fc — 1) of the columns C7tf> fc ~ 1-t- j) are independent. The 
columns C^' fc ~^ are all independent as a consequence of statement El among them, those 
with 1 < I < k — 1 can be expressed as linear combinations of CW> fe_1 ~ l ~-?) by eq. (jA.9j) . 
Only C^' fc ) and C^' -* are independent of all other columns because they are proportional 
to a degree-(fc — 1) polynomial in m, while all others are at most of degree (k — 2). Hence, 
only two independent vectors are added to the 2{k — l)-dimensional subspace spanned by 
CV> k ~ 1 ~ i ~i\ and the rank of is 

2(k - 1) + 2 = 2k. (A.22) 

It follows that each individual terms in the sum over m in eq. (j4.7j) depends only on 2k 
independent linear combinations of the coefficients gijk-i-j, 0, < i < k — 1, < j < k — i. 
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